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Abstract 

The formation of ZnO nanorings and nanohelixes with large polar surfaces observed in exper- 
iments [Nano Lett. 3, 1625 (2003); J. Am. Chem. Soc. 126, 6703 (2004)] is shown to be a 
result of the competition between elastic energy, spontaneous polarization-induced surface energy, 
volume energy, and defect-induced energy. It is found that nanorings and nanohelixes observed in 
experiments are stable and energetically favorable structures. 

PACS numbers: 62.25.+g, 77.65.-j 



1 



I. INTRODUCTION 



Piezoelectric material produces electrical charges when subjected to pressure, and 
stretches or contracts when an electrical field is applied,— which has been widely used as 
acoustic transducers in telephones, musical instruments, and so on. Piezoelectric crystals 
have no center of symmetry, which induces spontaneous polarization.— ZnO is a typical piezo- 
electric material with Wurtzite structure.^ Zn^+ and planes arranging alternately along 
its c-axis results in spontaneous polarization along this direction that induces divergent sur- 
face energy of Zn-(OOOl) and O-(OOOl) surfaces. To maintain them stable, some researchers 
propose that the surfaces are not reconstructed, in stead about one-fourth charges transfer 
between them.i^^ However, theoretical calculation and experiments reveal that the recon- 
struction is preferable to charge transfei^ in reducing the surface energy. 

After the discovery of ZnO nanobelts,— Kong and Wang synthesized novel ZnO nanobelts, 
especially, nanorings and nanohelixes whose surfaces were dominated by polarized Zn-(OOOl) 
and O-(OOOl) facets.— They attributed the observed phenomena as a result of surface po- 
lar charges and explained that the formation of these nanostructures as a consequence of 
minimizing the total energy contributed by spontaneous polarization and elasticity.— This 
model requires the ratio of the thickness t to the radius R of the nanorings to be smaller 
than a constant, and later more experimental data located in the green domain of Fig.Sa in 
Ref. Q supported this model. However, the optimal relation between the thickness and the 
radius of nanorings (black line of Fig.Sa in Ref. Illh in their concise model cannot fit well the 
experimental data. To improve this relation, an elaborate theory considering more factors 
than their model but remaining its main spirit is needed. 

Additionally, because of the symmetry of Wurtzite structure, we know that there are 5 
independent elastic constants for bulk ZnO that have been measured as cn = 190 GPa, 
Ci2 = 110 GPa, Ci3 = 90 GPa, C33 = 196 GPa, C44 = 39 GPa.— From these experimental 
values, we see Cn ~ C33, C12 ~ C13, and C44 ~ (cn — Ci2)/2. That is, bulk ZnO can be 
roughly regarded as an isotropic material, and thus its Young's modulus can be calculated 
as ~ 120 GPa from these experimental data. It is very surprising that Young's modulus of 
ZnO nanobelt is measured experimentally only ~ 50 GPa ,^^i^'^ which is much smaller than 
that of bulk ZnO. One needs to understand what makes this large difference of Young's 
modulus between bulk ZnO and the nanobelt. 
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In the present paper we try to address these issues. The rest of this paper is organized 
as follows. In Sec. |nj we construct the shape formation energy of ZnO nanobelts that 
consists of elastic energy, spontaneous polarization-induced surface energy, volume energy, 
and defect-induced energy. And then we put forward necessary conditions of structure 
existence. In Sec. IIIH we use necessary conditions of structure existence to explain nanorings 
and nanohelixes observed in experiments. We find that they are stable and energetically 
favorable structures. In Sec. llVt we give a brief summary and discussion. Some detailed 
derivations are supplied in Appendix. 

II. CONSTRUCTION OF SHAPE FORMATION ENERGY 

In order to quantitatively understand morphological problems of nanostructures in the- 
oretical point view, such as helixes of multi-walled carbon nanotubes^^*^ and cholesterol 
filaments,— Ou-Yang and his co-authors proposed a concept that the shape formation 
energy determines the morphology of structures. Their concept can be summarized and 
generalized as follows: a structure can exist if its shape formation energy is minimum and 
less than zero, provided that the thermal fluctuation can be neglected. We will define and 
construct the shape formation energy of ZnO nanobelts in this section, and then express 
explicitly the necessary conditions of structure existence. 

A. Definition of shape formation energy 

The synthesis process of ZnO nanobelts is briefiy described as follows . ^^'^^ When ZnO 
powders are heated to 1350 °C at a heating rate of 20 °C/min, they decompose into Zn^+ 
and O^^ at low pressure 10~^ torr. After a few minutes of evaporation and decomposition, 
the Ar carrier gas is introduced at a fiux of 25 standard cubic centimeters per minute. The 
condensation products are carried to a temperature zone of 400~500 °C and deposit onto 
an alumina substrate under the Ar pressure of 250 torr, where the end products — nanobelts 
are formed. This process can be simplified as a chemical "reaction" : 

BulkZnO ^ Condensation ^ ZnO Nanobelts. 

Schematic energy landscape of the synthesis process of ZnO nanobelts is shown in Fig. d 
where condensation products are regarded as a meso-phase. The shape formation energy F 
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is defined as the energy difference between ZnO nanobelt and condensation products with 
the same number of atoms. 

A ZnO nanobelt is a thin but long structure as schematically shown in Fig. |21 The 
symbols + and — represent Zn-terminated and 0-terminated surfaces, respectively. The 
curve C is its central line; ei is the tangent vector of curve C perpendicular to the cross 
section of the nanobelt; 62 and (parallel to [0001] direction of ZnO) are principal axes of 
the cross section. N is the normal vector of C that determines the bending direction of C in 
3-dimensional space, and 9 the angle between N and e.2- B is the binormal vector of curve 
C . w and t, the width and thickness of the nanobelt because of w ^ t in experiments,— >^ 
represent the dimensions in 62 and 63 directions, respectively. Because w and t are much 
smaller than the length of the nanobelt, the geometry of the nanobelt is determined uniquely 
by its central line C and the angle 6. Since the bending with 6 = ±tt/2 is preferred to others 
because of t <^ w, we concentrate on this case in this paper. The final form of the shape 
formation energy of the nanobelt should depend merely on geometric quantities of C, such 
as the curvature n and the torsion r. We propose that it consists of elastic energy, surface 
energy, volume energy, and defect-induced energy as discussed below. 



B. Elastic energy, surface energy, and volume energy 



The elastic energy per unit arc length can be written in an invariant quadratic form 
including infinitesimal displacements of the frame {01,62,63}: 



dei 

ds 



62 + 



k 



13 



ds 



■ 63 



+ /3 ( -3- ■ 63 

ds 



(1) 



where ds is the arc length element, {dei/ds) ■ 62 and (dei/ds) ■ 63 represent, respectively, the 
bending deformations around 63 and 62, while {de.2/ ds)-e^ represents the torsion deformation 
around 61. ki2 = Yotw^/12, ki^ = Yowt^ /12, (3ki^ with (3 = 2/{l + v) are the bending and 
torsion rigidities, where Yq and v are Young's modulus and poisson ratio of ZnO. For 
normal materials < z/ < 1/2, (3 is in the range between 4/3 and 2. 

Consider the geometric relations 62 = Ncos0 — Bsin^^, 63 = Nsin^ + Bcos^, and the 
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Frenet formula^ 



dei I ds 
dN/ds 
^ dQ/ds y 



K 
-K r 
-r 



N 



(2) 



Eq.(0) is simplified as 



k 



E,= -^[[a cos^ e + sin^ ^) fi:2 + /3 - d^/ds)^] 



(3) 



where a = kyilkx^- Especially, for Q = ±7r/2, the elastic energy per unit arc length can be 
written as: 

Ec = ^ + Pr') , (4) 

where kc = ki^ = Yowt^/ 12. 

Because the nanobelt has two large polar surfaces, we must take into account their surface 
energy. The surface energy per unit arc length is phenomenologically expressed as 



(5) 



where 7 is a constant quantity with the dimension of energy per unit area. 

The volume energy comes from two sources: One is the Gibbs free energy difference 
between the solid phase and the meso-phase (condensation products); another is the electric 
energy induced by the spontaneous polarization. The volume energy per unit arc length is 
phenomenologically expressed as 

Eg = gowt, (6) 
where go is a constant quantity with the dimension of energy per unit volume. 



C. Defect-induced energy 

As mentioned in Sec. HI the value of experimental Young's modulus of nanobelts (~50 
GPa) is less than half of bulk ZnO's value (~ 120 GPa). In order to compare the elastic 
constants of perfect ZnO thin films with bulk ZnO, we adopt a slab model and density 
functional theory (DFT) method to calculate them.^^ To reduce the computational task, 
we take the frozen-ion approximation and concentrate on the elastic constant cn. The 
computed ratio cf^^ /c^^^ as a function of the slab thickness is shown in Fig. El It is clear 
that the values of cf^^ /(^^^ are close to 1, which implies that the mechanical properties of 
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perfect ZnO nanobelts have little difference from bulk ZnO. Thus the nanobelts observed in 
experiments must contain defects that reduce Young's modulus of nanobelts significantly. 

There are several clues for the existence of defects in real nanobelts. A direct one is high 
resolution images of nanobelts. For example, the nanoring shown in Fig.4d of Ref. Q is not 
an perfect, uniform structure. High-resolution transmission electron microscopy also reveals 
that planar defects exist in the nanobelts . 

It is hard to write an exact expression of defect-induced energy. But if we only consider 
structures with \Kt\ <^ 1 and |rt| ^ 1, the defect-induced energy can be locally expanded 
up to the second order terms of nt and rt. It is expressed formally as 



where Ai, A2, and fj, are constant quantities with dimension energy per unit volume. To 
avoid the symmetry breaking of chirality, the linear term of rt is not included in the above 
expression. Comparing Eq.© with Eqs.(@]) and ®, we find that Ai and A2 play a similar 
role with the Young's modulus while A^ influences the value of volume energy density Qq. 
The rest term can induce a local spontaneous curvature of nanobelts randomly. Recently, 
the idea that defects induce a local spontaneous curvature has also been proposed to explain 
the cyclization phenomenon of short DNA.— 

D. Necessary conditions of structure existence 

The final form of the shape formation energy is thus written as 



where Y and g are the effective Young's modulus and the volume energy density of ZnO 
nanobelts, respectively. For example, the value of Y can be taken as the experimental value 
50 GPa. 

Now we estimate the effect of thermal fluctuation on the shapes of nanobelts. Considering 
the typical growth temperature T ~ 10^ K, the Young's modulus Y = 50 GPa, the typical 
thickness t ~ 10 nm, and the width w > 50 nm in the experiments,— »ii we estimate the 
persistence length of nanobelts as Ip = Ywt^ /12T ~ 1 cm,— which is much larger than the 



Ed = - [AiiKtf + A2(rt)2 + fx{Kt) + A3] wt 



(7) 



F 
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typical length of nanobelts of several microns. Thus the effect of thermal fluctuation can be 
neglected and the necessary conditions of structure existence can be expressed as: 



5F = 0, 
5^F > 0, 



F < 0, 



(9) 
(10) 

(11) 



where 6F and represent the first and second order variations of the shape formation 
energy, respectively. If a structure satisfies the above three conditions, we call it stable and 
energetically favorable structure. 

From Eq.(jni), we arrive at the following shape equations of nanobelts: 



KT -\ 

2 



4r 



K 



+ ^ + 3kt^ 

K 



+ 



t 



XT 



0, 



2t3 



0, 



(12) 
(13) 



where x = 12/i/F, t] = 487/F, and ^ = 2Ag/Y. ()s represents the derivative respect to 
s. The detailed derivations of above two equations are shown in Appendix. The necessary 
conditions of structure existence are transformed into Eqs. (|Tn|) - (fT!?jl . 



III. EXPLANATION TO EXPERIMENTS 



Now let us use Eqs. fll()|l - pH|l to explain nanorings and nanohelixes observed in experi- 
ments. 



A. Nanorings 

For nanorings with radius one has k = 1/R and r = 0. Eq. ()12j) is transformed into 

it/Rf = at/R) + v/R, (14) 

while Eq. (llH|l is trivial. By fitting this equation to the experimental data^i as shown in 
Fig. H we obtain ^ = -0.001 ± 0.0007 and rj = 0.03 ± 0.01 nm. 

Now we test condition ()10|) . In polar coordinate system {p,4>), any perturbation in the 
vicinity of a cycle with radius R can be expressed as p = R{1 + Xln^oo '^n^*""^)' where 
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a-n = a* (n = 0, 1, 2, ■ ■ ■ , oo) are small constants. The second order variation is evaluated 
explicitly as 

S^F = ^(n^ - iy\an\^Ywtyi2R > (15) 

n 

under the validity of Eq. (fT^ . 

Next, nanorings observed in experiments must satisfy condition (jllll . This requires F = 
T{Ywt^(t/ R—x)/^ < to be valid for all experimental data. Thus we can take the supremum 
of t/R in Fig. HJas the value of x, i-e., x = 0.08. 



B. Nanohelixes 



For a nanohelix with radius tq and pitch p, the pitch angle is = arctan(j9/27rro). The 
curvature and the torsion are k = cos^ip/vQ and r = — simp cos if/ tq, respectively. Eq. p2|) 
is transformed into 

cos V9 + (3p — 2) sm ipcos ip -\ sm (p = — — 5 , (16) 

while Eq. fjl3p is trivial again. Eq. ()16p sets a relation between the thickness, the pitch angle 
and the radius of nanohelixes, from which we can estimate the thickness. The results by 
using the data in Ref. Q are summarized in Table HI These values are close to the typical 
value 10 nm observed in experiments.— 

Additionally, the shape formation energy can be calculated as 

„ Ywt^L 



12r2 



COS 09 H sm 2o9 cos 2o9 

4 t 



(17) 



where L is the total length of the nanohelixes. Its value is listed in Table d for three 
nanohelixes observed in the experiments using x = 0.08. Moreover, we also found the 
matrix 

d'^F/drl d'^F/dr^d^ \ 
d^F/drod^ d^F/dip^ j 
is positive definite for these nanohelixes. Therefore, these nanohelixes are stable and ener- 
getically favorable structures. 
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IV. CONCLUSION 



In the present work, we construct the shape formation energy of ZnO nanobelts that con- 
sists of elastic energy, spontaneous polarization-induced surface energy, volume energy, and 
defect-induced energy. We put forward the necessary conditions of structure existence and 
find that nanorings and nanohelixes of ZnO observed in experiments are stable and energet- 
ically favorable structures. We obtain an important result that nanorings and nanohelixes 
must satisfy Eqs.lfT^ and (fT^. respectively. We notice that it is easy to measure experi- 
mentally the values of the radius of nanorings or the pitch and radius of nanohelixes, while 
it is difficult to measure their thickness of about ten nanometers in high precision because 
the error of InLens detector is ±2 nm .^°i^^ Eqs. (ll4|) and provide a way to overcome 
this difficulty. The present theory may also be developed to explain the recent synthesized 
sup erlat t ice- st r uct ured nanohelxes .— 

Because of the influence of defects. Young's modulus Y and volume energy density g are 
quite different from the value of perfect and uniform nanobelts. For example, Y = 50 GPa 
is just half of the value of perfect nanobelts. g is estimated as ~ —10^ J/m^ by ^ = —0.001. 
This value departs largely from the cohesive energy ~ —10^° J/m^ of ZnO.— The surface 
energy density 7 is estimated as 0.03 J/m^, which cannot compare with the cleavage energy 
~ 4 J/m^ for perfect Zn-(OOOl) and O-(OOOl) surfaces obtained by ab initio calculations.— 
There is no theory on the parameter ^ at present, which is expected to depend on the 



defect density and experimenta 



conditions. We estimate its value fi ^ 10^ J/m^ from the 



experiment published in Ref . ^ . 

In the above discussion, we do not include flexoelectric effec t ^^1^^ of piezoelectric materials 
in nanoscale. This effect is too small for crystals^— to influence the above results in statics. 
But it may play a role in the kinetic growth of ZnO nanobelt so that nanorings with Zn- 
terminated and 0-terminated inner surfaces have different weights in the end products. 
Through analyzing the weights, one may reveal the flexoelectric effect of ZnO in experiments. 
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APPENDIX 

We present the derivation of Eqs. (jl2|) and (jl3|) from the first order variation of functional 
(jEJ under the periodic boundary or fixed boundary. The basic idea is similar to Ref. Isfll 
When we perform the variational calculation, the volume and density of a ZnO nanobelt 
are allowed to vary, while the total mass is conserved. However, the variation is performed 
for infinitely small perturbation which has very small effect on the elastic coefficients in 
Eq.®. The final Euler-Lagrange equations obtained from the variational process including 
this effect are the same as those derived from the variation with elastic coefficients fixed to 
constant. 

Let El = Bi, E2 = 'N , = B and dEi = WijEj, where Einstein's summation convention 
is used and "c?" is an exterior differential operator. Thus we have dr = Eids, W12 = 
hids, Wis = 0, W23 = rds by considering the Frenet formula. 

First, we consider 6r = Q2E2 and 6Ei = QijEj. Using the formula d6r = 6dr and 
d6Ei = 5dEip^ we have 

5ds = 1^2^21 = —i^dsfl2, (18) 
^12^5 = dn2, (19) 
^^13 = rn2, (20) 
5wij = dVLij + VLikWkj - Wik^kj- (21) 

Because 6wi3 = 0, we obtain 

n23ds = [2d (rfia) - f^2C^r]/K. (22) 

Using above five equations, we can prove 

6Kds = dni2 + {k^ - T^) VL2ds, (23) 
^Tds = dil23 + 2KTdsVL2. (24) 

Considering Stokes theorem and above seven equations, we arrive at 

5 j K^ds = y 2 {f^ss + /^V2 - f^r^) ^2ds, (25) 
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6 J T^ds 



At 



2tI 



K J s K 



(26) 
(27) 



Consequently, we obtain 

Ywt^ 



SF 



K 

y 



KT 



+ 



12 

wlfifr'^ - (27 + gt)]KVL2ds. 



4 



4r 



2r2 



(28) 



Because VL2 is an arbitrary function, 5F = gives Eq. ()12j) . 

Secondly, we consider 6r = fl-sE^. Using the formula dSr = Sdr and d6Ei = 6dEi, we 
have 



Sds 

flisds 
5wii 



ilswsi = 0, 
dfl3, 



kj- 



Because Swis = 0, we obtain 



^23ds = [d^li3 — r'^risds]/ K. 



Using above five equations, we can prove 

Suds = —2d (tQ^) + Q^dr, 
drds = dQ23 + ndQ^. 

Considering Stokes theorem and above seven equations, we arrive at 



5 J K^ds = J {AksT + 2kTs) n^ds, 



6 J T^ds 



2 

TgQsds. 



T n - ( TK ) - — 

^ K/ss 



Consequently, we obtain 
Ywt^ 



5F 



12 



2KsT + KTs + P 



T^T - [TK] 



24fiTs 
Yt 



Q'^ds. 



(29) 
(30) 
(31) 
(32) 



(33) 



(34) 
(35) 



(36) 
(37) 
(38) 

(39) 
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Because is an arbitrary function, SF = gives Eq.(fT^. 
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TABLE I: Radii, pitches, estimated thicknesses and shape energies for different hehxes. The data 
for radii and pitches are taken form Ref. lid . L represents the length of hehxes. 



ro (nmj 



p (nm) 



t (nm) 



F {Ywt^L/l2rl) 



342 
175 
240 



222 
133 
380 



12.4~12.5 
8.5~8.6 
9.6~9.8 



-1.1 
-0.6 
-0.7 



Energy 




BulkZnO 



FIG. 1: Schematic energy landscape of the synthesis process of ZnO nanobelts. 




t/2 



FIG. 2: Schematic figure of a ZnO nanobelt. C represents its central line, ei is the tangent vector 
of curve C that is perpendicular to the cross section of the nanobelt. e2 and 63 (parallel to [0001] 
direction of ZnO) are the principal axes of the cross section. N and B are the normal and binormal 
vectors of C, respectively. 9 is the angle between N and 62. 
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FIG. 3: Relative value of elastic constant cn for a slab and bulk ZnO. 
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FIG. 4: Relation between l/R and the ratio t/R of nanorings with radius smaller than 4000 nm. 



The data are taken from Fig. 5a of Ref. 
The dash line is the supremum of t/R. 
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The solid line is a fitting curve of the data by Ea. H14|) . 
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